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The accumulation of 7’Rn and ”°Rn progeny in poorly ventilated environments poses the risk of natural 
radiation exposure to the public. A previous study indicated that satisfactory results in determining the *Rn 
and *°Rn progeny concentrations by measuring the total alpha counts at five time intervals within 560 min 
should be expected only in the case of high progeny concentrations in air. To complete the measurement within 
a relatively short period and adapt it for simultaneous measurements at comparatively lower 7Rn and Rn 
progeny concentrations, a novel mathematical model was proposed based on the radioactive decay law. This 
model employs a nonlinear fitting method to distinguish nuclides with overlapping spectra by utilizing the alpha 
particle counts of nonoverlapping spectra within consecutive measurement cycles to obtain the concentrations of 
?2Rn and Rn progeny in air. Several verification experiments were conducted using an alpha spectrometer. 
The experimental results demonstrate that the concentrations of 7’Rn and 7°Rn progeny calculated by the 
new method align more closely with the actual circumstances than those calculated by the total count method, 
and their relative uncertainties are all within + 16%. Furthermore, the measurement time was reduced to 
90 min, representing an acceleration of 84%. The improved capability of the new method in distinguishing 
alpha particles with similar energies emitted from ”!ŝPo and *’’Bi, both approximately 6 MeV, contributed to 
realizing more accurate results. The proposed method has the potential advantage of measuring relatively low 
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concentrations of ?’Rn and 7°Rn progeny in air more quickly via air filtration. 
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I. INTRODUCTION 


Radioactive gases 222Rn and Rn are colorless, tasteless, 
and odorless. They are universally present in the walls, floors, 
and ceilings of buildings, as well as in the surrounding soil, 
before permeating into the cracks of a building and entering 
indoor spaces [1, 2]. Their decay products typically accumu- 
late in closed spaces or poorly ventilated areas such as base- 
ments, mines, and warehouses [3, 4]. The generated progeny 
atoms recombine with air ions and other air impurity clusters 
to form clusters of size 0.5 nm — 5nm, which are referred to 
as radioactive unattached particles [5—7]. After this process, 
the clusters attach to submicron-sized aerosol particles in the 
air within 1-100 s, forming radioactive aerosols or attached 
particles with sizes ranging from 100 nm to 500 nm [8, 9]. 
During breathing, radioactive progeny aerosols are inhaled 
into the respiratory organs. Some aerosols remain adsorbed 
in the respiratory system, where they continue to undergo 
radioactive decay [10]. During this process, the progeny of 
?22Rn and °Rn can cause continuous tissue damage, leading 
to lung cancer. The inhalation doses are contributed by their 
progeny, but not by the gas itself [11]. Therefore, numer- 
ous researchers conducted studies on detectors for measuring 
22Rn and ?®Rn progeny concentrations. 
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Various detectors have been developed for measuring al- 
pha particles emitted from 7”*Rn and 7”°Rn progeny in air, in- 
cluding ZnS(Ag)-coated scintillation cells [12, 13], the imag- 
ing plate (IP) [14], solid-state nuclear track detectors (SS- 
NTD) [15, 16], silicon semiconductor detectors [17, 18], and 
others. An alpha spectrometer, equipped with ULTRA ion- 
implanted silicon semiconductor detectors, features a high- 
energy resolution, low background, and appropriate detection 
efficiency. Therefore, it is frequently employed for the identi- 
fication and quantitative analysis of mixtures of radionuclides 
based on the alpha particles emitted during decay [19, 20]. 


Previously, methods based on air filtration followed by 
measurements using an alpha spectrometer have been widely 
explored for the determination of mixed ??Rn and ??°Rn pro- 
genies in environmental air [21—25]. Harley et al. proposed a 
method that includes a 60-min sampling period, with counts 
considered at five intervals within 2440 min, to estimate the 
individual concentrations of 7??Rn and 7°Rn progeny. Hence, 
the lung dose due to 7°Rn progeny [21] is assessed. A 
comparable method employed in high-background areas in 
Yangjiang County, China, involved 30-min sampling, fol- 
lowed by five measurement intervals within 560 min, and 
ensuring on-site measurement work completion within ten 
hours [22]. For more accurate measurement of the concen- 
trations of mixed 7?’Rn and 7°Rn progeny, a set of optimum 
time intervals was obtained by comparing the uncertainties in 
the total alpha spectroscopy data across the three sets of mea- 
surement intervals. Each set of time intervals lasted 560 min 
and was calculated using the weighted least-squares method 
[23]. To assess the reliability of this air-filtration method, the 
sensitivity of the measurements of 7”’Rn and 7?°Rn concen- 
trations in air to variations in alpha counting at three and five 
intervals (560 min) was examined [24]. They indicated that 


the ?”’Rn progeny reached saturation activities on the filter 
after 3 h of filtering, whereas the saturation activities of 7?°Rn 
were achieved after three days owing to the different half- 
lives in the °Rn chain. The influences of ???Rn and Rn 
concentrations, filtering duration, and choice of measuring 
interval (560 min) on the relative standard deviations were 
further analyzed [25]. Their study indicated that satisfactory 
results using the total count method should be expected only 
in the case of high radon and thoron progeny concentrations 
in air. 

In contrast to the total alpha-counting method, which re- 
quires five counting intervals, the alpha spectroscopy method 
enables the direct differentiation of alpha particles with vary- 
ing energies emitted by 7!*Po, 7!4Po, *!?Bi, and 7!*Po. This 
is widely preferred, particularly in the context of modern re- 
quirements for high measurement accuracy [26-29]. In 1978, 
Kerr introduced an a spectrometry technique to measure the 
activity concentration of short-lived radon and thoron proge- 
nies. This method, termed as Kerr—Tn, utilizes an inverse ma- 
trix and obtains five alpha counts within the region of interest 
(ROD. Three counting intervals were established, each with a 
sampling time of 10 min: 2—12 min, 15-30 min, and 200-220 
min. The intervals were initiated after the completion of sam- 
ple collection [30]. Compared to the Kerr-Tn method, both 
the PKU-Tn and Wicke-Tn methods exhibit the advantage 
of synchronizing the first counting interval with the sampling 
time. With an extended 30-min sampling time and a total 
measurement duration of 180 min, this results in an increase 
in methodological sensitivity by nearly 8.97 times and 3.56 
times, respectively [31]. However, adopting a method of si- 
multaneous sampling and measurement may lead to broad- 
ening of the measured alpha energy spectrum because there 
is a need for a distance between the detector and filter mem- 
brane to allow for gas circulation. This in turn affects the 
accuracy of the measurements [32]. Furthermore, previous 
fitting methods demonstrated significant advantages in alpha 
energy spectrum measurements. For the deconvolution of the 
alpha-particle spectra, improved peak shape formulae were 
implemented as fit functions in a spreadsheet application, and 
the optimum fit parameters were searched using built-in op- 
timization routines [33]. Considering the entire spectrum 
background, background noise, and peaks with low statisti- 
cal counts, a model was developed by fitting the typical shape 
of the alpha peak [34]. 


To complete the measurement within a relatively short pe- 
riod of time and adapt to simultaneous measurements at com- 
paratively lower *”’Rn and 7°Rn progeny concentrations, 
a new mathematical model based on the radioactive decay 
law was proposed by performing a non-linear fitting on the 
counts of alpha particles with different energies over contin- 
uous measurement periods to deduce the concentrations of 
?22Rn and >°Rn progeny in air. Several measurements were 
performed using the alpha spectrometer. The results of the 
proposed method are then compared with those of the total 
count method [25]. This indicates that the newly proposed 
method shows more potential advantages in the fast simulta- 
neous measurement of relatively low concentrations of 7’Rn 
and 7?°Rn progeny in air. 


Il. MATERIALS AND METHODS 
A. alpha emitter from Rn and °Rn progeny 


The decay chains of ”’Rn and 7°Rn are shown in Fig. 
1 [35, 36], excluding branches with a low probability. There 
are some short-lived alpha emitters of radon and thoron 
progeny in the air, including RaA (7!*Po, 6.0 MeV), RaC’ 
(Po, 7.69MeV), ThA (7!Po, 6.78 MeV), ThC (7!°Bi, 
6.05 MeV), and ThC’ (7!*Po, 8.78 MeV). Given the air filtra- 
tion method employed in this study, progeny with extremely 
short half-lives, such as 7!*Po and *!*Po, quickly reach equi- 
librium with their parents during the sampling process. Ad- 
ditionally, it is assumed that 7!?Pb originates directly from 
20Rn due to the short half-life of 7!°Po. In reality, °! Bi 
emits an alpha only 36% of the time, whereas the remain- 
ing 64% of the time it emits a beta. This is followed almost 
immediately by an alpha through the decay of 7!*Po (with a 
half-life of 0.3s). Furthermore, 7!*Po decays with an alpha 
energy of 8.78 MeV and 7°8T] decays via beta decay only. It 
is impossible to distinguish the activity on the filter because 
of 6.05 MeV alpha particles from ThC and 6.0 MeV alpha 
particles from RaA. However, the activity can be partitioned 
between RaA and ThC using the 8.78 MeV alpha particle ac- 
tivity from ThC’ [30]. 


poo e eee ea a a 


ern d as $7 eo: @- ne 
(a 5.49Mel 0. 6.00Mev ia 


ThC’ ` 
ThC pon eA a 1 
55.68 s 10.64 h 4 
Pep en i 
.29MeV a S gid 1 
1 ma ; 
7 1 
t 1 
% 
S ThC” ot 


Fig. 1. (Color online) Decay chains of *’Rn and Rn 


B. Measurement method of the progeny of 7”Rn and ”°Rn 
via alpha spectroscopy 


The new method of simultaneously measuring the progeny 
concentrations of ??Rn and ”°Rn by using an alpha spec- 
trometer includes three time intervals: sampling duration 
(Tı), time interval (T2) between the end of sampling and be- 
ginning of the measuring the progeny of 7”’Rn and 7°Rn 
on the filter membrane, and time duration (T3) for alpha 
spectrometer measurements from the beginning to the end. 
Three assumptions are made: 1) the flow rate remains steady 
during the sampling process; 2) the concentrations of RaA, 
RaB, RaC, ThB, and ThC in the environment remain constant 
throughout the sampling process; and 3) the filter exhibits 
the same collection efficiency for different °? Rn / 220Rn pro- 
genies, and its self-absorption can be considered negligible. 
This procedure is illustrated in Fig. 2 and can be described 


as follows: 
(1) Sampling 


Before the sampling begins, the number of atoms in ??Rn 
and 7?°Rn progenies on the filter membrane is zero. The 
switch of the vacuum pump is turned on to draw air contain- 
ing the progeny of 7”’Rn and 7”°Rn into the sampler through 
a flow meter, where 2”*Rn and ?®Rn are collected on the fil- 


dNraa/dt = QGCRaa/ARaA 


ter. Subsequently, air is returned to the chamber with mixed 
concentrations of °”Rn and ”°Rn. The sampling duration 
lasts for fifteen minutes. During the sampling process, the 
quantity of *??Rn and ?°Rn progeny on the filter membrane 
continuously increases. However, this quantity also decreases 
owing to the decay of parent nuclides and their decay. This 
can be expressed using the following system of differential 
equations. 


— NradARad 


dNpap/dt = OGCraB/ARaB + NRaA ÀRaa — NraBrARaB 


dNrac/dt = QGCRac/ARac 
dNrnp/dt = QGCrnp/AtThB 
dNrnc/dt = QGCrnc/Atac 


where G is the filtration efficiency, Q is the flow rate, CRaa, 
Crap, Crac, Crnp, Crnac are the activity concentrations 
of RaA,RaB,RaC,ThB and ThC, respectively. ARaa, ARaB, 
ARaC; ATHB; Arac are decay constants. Nraa, NRaB, Nrac, 
Nrnp, Nrnc are the number of atoms of RaA,RaB,RaC,ThB 
and ThC on the filter membrane during the sampling process. 

The solution to Eq. | shows the relationship between the 
accumulated number of atoms of RaA,RaB,RaC,ThB and 
ThC on the filter membrane at the end of the sampling time 
T and the activity concentrations of RaA, RaB, RaC, ThB, 
and ThC in air. 

(2) Filter membrane processing 

After completing the sampling, the filter membrane was 


+ NRaBÀRaB — NracArac (1) 
— NThBÀTHB 
+ NThBAThB — NrncAtnc 


removed and placed in the sample cell of the alpha spectrom- 
eter for measurement with a one-minute time interval. In this 
process, the 20Rn gas did not remain on the filter membrane, 
and its progeny 7!°Po on the filter membrane completely de- 
cayed. Therefore, there was no spectral interference from the 
alpha particles emitted by the decay of 7’Rn and 7!°Po dur- 
ing the measurement. At the end of sampling, the quantity 
of progeny of 7”Rn and Rn on the filter membrane be- 
gan to decrease because of decay. Meanwhile, the quantity 
of progeny nuclides also increases owing to the decay of the 
parent nuclides. This can be expressed by the following dif- 
ferential equation: 


dNraai/dt = —NRaai ARaA 

dNpapi/dt = NRaA1ÀRaa — NRaB1ÀRaB 

dNraci/dt = NRaB1 ÀRaB — NraciArac (2) 
dNrnpi/dt = —NrnpiAtns 

dNrnci/dt = NrnpidAtns — NrnciAtac 


where Nraat, Nrapi, Nraci, Nrnpi, and Nrnci denote 
the number of atoms of RaA, RaB, RaC, ThB, and ThC, re- 
spectively, on the filter membrane during time interval T>. 

The solution to Eq. 2 establishes a relationship between the 
number of atoms of RaA, RaB, RaC, ThB, and ThC on the 
filter membrane at time Tə and the number of atoms of RaA, 
RaB, RaC, ThB, and ThC at the end of the sampling period. 

(3) Measurement process 

After processing, the filter membrane was sent to the in- 
ternal cell of an alpha spectrometer. The alpha spectrome- 
ter was connected to a computer via a data cable, and the 
measured counts of different alpha particles under a vacuum 
pressure below 1000 mTorr were displayed in real time using 
MAESTRO software on the computer. At the beginning of 
alpha spectrometer measurement on the filter membrane, the 


number of atoms of RaA, RaB, RaC, ThB, and ThC on the 
filter membrane are represented as NRaaz2, NRaB2, NRac2, 
Nrnpe2, and Nraco, respectively. According to Eq. 1 and 2, 
the activity concentrations of RaA, RaB, RaC, ThB, and ThC 
in air can be deduced. After the alpha spectrometer mea- 
surement was initiated, the progeny of 7”’Rn and 7°Rn on 
the filter membrane continues to decay according to Eqs. 2. 
During the measurement process, the number of RaA, RaB, 
RaC, ThB, and ThC atoms on the filter membrane at any given 
time t are represented as NrRaa3, NRaB3, NRac3, Nrnp3, and 
N-rnc3, respectively, and they can be determined as follows: 


Nraaa(t) = Nraaze *R™4* (3) 
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Fig. 2. (Color online) Process of simultaneously measuring the progeny concentrations of *’Rn and *°Rn using the alpha spectrometer 
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During the alpha spectrometer measurements, the activities 
of RaA, RaB, RaC, ThB, and ThC on the filter membrane are 
denoted as Agaa, ÁRaB, ARac, ATnp, and Arno, respec- 
tively. Distinguishing between the activity on the filter at- 
tributed to 6.05 MeV alpha particles from ThC and 6.0 MeV 
alpha particles from RaA is unfeasible. However, the activ- 
ity can be apportioned between RaA and ThC by utilizing 
8.78 MeV alpha particle emissions originating from ThC’. 
According to the radioactive decay characteristics of 7’Rn 
and ??°Rn, the alpha particle counts of 8.78 MeV characteris- 
tic peak originate from 0.64 x Arnc; the alpha particle counts 
of the 7.69 MeV characteristic peak originate from Arac; 
and the alpha particle counts of the 6.0 MeV characteristic 


NThC (i) _ . 1 
T = 0.64E Atno (ir z7) 


( NThB2ÀThB 


N: grae (12) 
ÀThC — ÀThB mcz) 


peak originate from Apa, + 0.36 x Arc- 

Short measurement cycles were implemented with a mea- 
surement period T, and there are ¿ measurement cycles. The 
total measurement time is T3. In the i-th measurement 
cycle, the net alpha counts are obtained in their respec- 
tive alpha spectrum of regions of interest, denoted by ROI- 
1 (8.78 MeV), ROI-2 (7.69 MeV), and ROI-3 (6.0 MeV), 
and they can be represented by nracli), NRac(i), and 
NRaa+ThC (i). Owing to the short duration of each measure- 
ment cycle, the average counts for the ith measurement cycle 
are approximately equal to the values of Apnc, Arac, and 
ARaa at the midpoint of the measurement cycle, and they can 
be expressed as: 
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where F denotes the detection efficiency of the detector for 
alpha particles, i denotes the counting cycle of alpha particles, 
and iT — 1/2T denotes the midpoint time of the i-th counting 
cycle. 


Based on Eqs. 13, 14, and 15, the number of atoms 
(Nraa2, Nrapo, Nrace, Nrhg2, and Nyc) can be deter- 
mined through the nonlinear fit of the experimental data, and 


( NThB2ÀThB 
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Athc — ATHB moz) i 


(5) 


the calculation processes are as follows: 


(1) NThg2 and NrThco can be obtained using Eqs. 13 
to perform a nonlinear fitting of the experimental data for 
(nrnc(i))/t. 

(2) By substituting the value of Nrhc2 into Eqs. 15 
and employing nonlinear fitting of the experimental data for 
(ngaa (i))/t using Eq. 15, Nraag can be obtained as: 


(3) By inserting Npaag into Eq. 14 and utilizing the non- 
linear fitting of the experimental data for (NRaA+Tnc)/t us- 
ing Eqs. 14, NRag2 and Nracg are determined. 

After obtaining the values of Nraao, Nrapo, Nraco, 
Nrnpo2, and Nyncg, the activity concentrations of RaA, RaB, 
RaC, ThB, and ThC in air can be deduced by reversing Eqs. 
l and 2. 


C. Limit of detection and error analysis 


The net signal level (instrument response) above which an 
observed signal can be reliably recognized as "detected is 
defined as the critical limit (Lc). This indicates that a par- 
ticular nuclide is definitely present in the sample [35]. The 
same values are defined in ISO standard 11,929, but the math- 
ematical assessment of ionizing radiation measurement uti- 


where oc, (i=1,2,3,4,5) denotes the uncertainties of activity 
concentration of RaA, RaB, RaC, ThB, and ThC, og denotes 
the standard uncertainty of the detection efficiency, oq de- 
notes the standard uncertainty of the filtration efficiency, and 
aQ denotes the standard uncertainty of the flow rate. 


I. RESULTS AND DISCUSSION 


Several experiments based on air filtration of Rn and 
?20Rn progeny concentrations were conducted in mixed 7”’Rn 
and ??°Rn chamber at the University of South China. This 
chamber is made of stainless steel and has an inner effective 
volume of 3m. It is connected to the gas path, allow- 
ing continuous compensation for 7??Rn and 7?°Rn progeny. 
The flow rate for each sampling experiment was calibrated 
to 15 + 0.45 L- min“! using a soap bubble flowmeter, and 
the sampling duration lasts for fifteen minutes. The volume 
of the sampled gas occupies only 8% of the entire chamber 
volume, and continuous progeny replenishment ensures rela- 
tive stability in the measurement environment. The sampling 
filter membrane is a 0.8 wm Millipore AA-type microfiltra- 
tion membrane with a diameter of 25 mm. The filtration effi- 
ciency of the filter membrane is 0.990 + 0.002, and its self- 
absorption is considered negligible. Aerosols are generated 
using an aerosol generator, and an aerosol electrometer, pro- 
duced by TSI Incorporated, is used to monitor the number 
concentration of aerosols. The aerosol concentration varied 
from 8 x 103 cm~? to 10 x 10° cm~? in the experiments. 

The measurements were performed at 27 degreeCelsius 
and 75% relative humidity (RH) using an alpha spectrome- 
ter. This alpha spectrometer, manufactured by ORTEC Cor- 
poration, utilizes an ion-implanted surface-passivated silicon 


lizes Bayesian statistics, which comprise calculations with 
conditional probabilities [37-40]. According to ISO standard 
11,929, the decision threshold or critical level (Lc), assum- 
ing small number of background counts (C},), was calculated 
using the following formula: 


Lo( in counts ) = k - y2 - (Cy + 1) (16) 


where k denotes the confidence factor (k=1.645 at confidence 
level of 95%). 


The uncertainties in the activity concentration(ac, ) can be 
expressed by Eq. 17 [37-40]. It is composed of the un- 
certainty of counts Ni, detection efficiency E, filtration effi- 
ciency G, and flow rate Q, whereas the uncertainties of sam- 
pling and counting times are usually ignored. 
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detector with an 900 mm? probe. It realizes an energy resolu- 
tion of 20 keV under a vacuum pressure below 1000 mTorr. 
Under these measurement conditions, the alpha particles ex- 
hibit a longer range and minimal energy loss with a measure- 
ment distance of only 1 mm and no spectral broadening oc- 
curs. Consequently, the distinct characteristic energy peaks of 
222Rn and Rn progenies can be clearly distinguished using 
1024 channels. The alpha spectrometer calibrates its alpha 
detection efficiency (E) using a *4'Am electroplated surface 
source. The resulting mean value for FE was 33.3%, accom- 
panied by 2.0% uncertainty. 

To compare the accuracy of the results obtained from the 
proposed model with those of the total count method pre- 
sented in [25], the measurement cycle of the alpha spectrom- 
eter was set to 1 min and lasted for more than 560 min. This 
method can provide alpha counts at different time intervals 
as required by both methods. The automatic measurement 
process using an alpha spectrometer was realized by writ- 
ing JOB files in MAESTRO software, which can display and 
store real-time measurement data. A portion of the measured 
data is presented in Fig. 3. 

Considering the impact of statistical fluctuations on net al- 
pha counts and the influence of measurement intervals and cy- 
cle periods on the standard error of the calculated results, nine 
short counting cycles (i) were selected consecutively from the 
beginning of the measurement, with each counting cycle (T) 
set to 10 min. The background count for a measurement time 
of 90 min was 31 + 3. According to Eqs. 17, the critical 
limit (Lc) of the alpha spectrometer for the counts obtained 
using the new method was 13. The net counts of alpha par- 
ticles emitted from ?”Rn and ”°Rn progeny were recorded 
based on alpha particle energies of 8.78 MeV, 7.69 MeV and 
6.0 MeV, and their ROI-1 (8.78 MeV), ROI-2 (7.69 MeV), 
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Fig. 3. (Color online) Alpha counts with energies of 6.0 MeV, 7.69 MeV and 8.78 MeV measured via display in Maestro software. 


and ROI-3 (6.0 MeV) were 2-6.3 MeV, 6.3-8 MeV, and 8- 
10 MeV, respectively. 


Origin software (OriginLab Corporation) enables users to 
customize functions for nonlinear fitting, ensuring the accu- 
rate modeling of complex data relationships. This capabil- 
ity provides enhanced precision and adaptability in nonlinear 
data analysis, making it a widely used tool for calculating 
22Rn and 7°Rn concentrations [41, 42]. The fitting method 
is based on the radioactive decay law that uses nonlinear fit- 
ting by Origin2021 and return the best-fit parameters and their 
statistical errors for each free parameter fitted. In this study, 
the fitted parameters NRaa2, NRaB2, NRac2, Nrnpe2, and 
Nrnce are the primary focus. Fig. 4 illustrates three experi- 
mental results for the net alpha particle counts measured by 
the alpha spectrometer, which varies with the measurement 
period. The alpha counts from RaA and ThC, both around 
6 MeV, were counted, and the results of the three experiments 
are shown in Figs. 4b, e, and h. Given significantly shorter 
half-life of RaA when compared with that of ThC, their alpha 
counts initially decreased rapidly, followed by a slow rise. 
The alpha counts from ThC and RaC are presented in Figs. 
4a, (d), and (g), and Figs. 4c—(i), respectively. According 
to Eqs. 13, 14 and 15, the number of atoms of RaA, RaB, 
RaC, ThB, and ThC on the filter membrane (Nraa2, NRaB2, 
NRac2, Nrnp2, and Nypce) at the initiation of alpha spec- 
trometer measurement was derived through non-linear fitting 
of alpha counts with energies of 6.0 MeV, 7.69 MeV, and 
8.78 MeV using Origin2021 software. By inserting the values 
of NRaa2, NRaB2, Nracz, Nrnp2, and Nraco into Eqs. 1 
and 2 using MATLAB, the concentrations of 7”’Rn and 7”°Rn 
progenies in the air can be obtained, and the detailed data are 
listed in Table 1. 


Only minor changes were observed in the number of 
aerosol particles and environmental parameters, and the 
sources of 7*Rn and Rn remained unchanged during the 


Table 1. Three experimental results of simultaneous measurements 
of 7’Rn and Rn progeny concentrations (k=1). 


the number of 


i doa Progeny atoms onthe R? Activity snarl 
pil filter membrane mea aa) 
1 Nraa2 229+ 19 0.93 17.234+1.537 
Nrag2 2270+ 133 0.96 3.845+0.364 
Nrac2 1220+97 0.96 2.809+0.362 
Nrnp2 30955 + 2244 0.67 2.539+0.199 
Nrnc2 1646+ 123 0.67 1.298-+£0.124 
2 Nraa2 308+ 46 0.63 23.175+3.495 
Nrag2 40534124 0.99 7.569+0.491 
Nrac2 2499+93 0.99 6.098+0.384 
Nrnp2 60926 + 4340 0.80 4.997+0.386 
Nrnc2 23654229 0.80 1.729+0.225 
3 Nraa2 275445 0.65 20.665+3.453 
Nrag2 24704122 0.97 4.006+0.446 
NrRac2 1358+89 0.97 3.178+0.338 
NrThge2 38466 + 1593 0.92 3.155-£0.161 
Nrne2 1550+ 84 0.92 1.146+0.088 


process of conducting the three measurement experiments. 
Table 1 shows that R? for Nrhgz and Nrpco in the first ex- 
periment, as well as for Nraag in the second and third exper- 
iments, are not satisfactory. This may be due to the impact 
of statistical fluctuations in counts over consecutive cycles. 
However, based on the results of the three experiments, the 
differences observed in the concentrations of the same types 
of ?”’Rn and Rn progeny calculated using the same method 
in the three experiments may be attributed to minor variations 
in the quantity of aerosol particles suspended in the environ- 
ment and other environmental parameters. The poor good- 
ness of fit of one of the nuclides slightly impacts the overall 
experimental results. 


The 7”*Rn and Rn progeny concentrations calculated us- 
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Fig. 4. (Color online) Three experimental results of net alpha particle counts measured by the alpha spectrometer (with energies of 8.78 MeV, 
6.0 MeV, and 7.69 MeV) vary with the measurement period: the first experiment (a, b, c), second experiment (d, e, f), third experiment (g, 
h, i). i represents the i-th measurement period, with each cycle lasting 10 min. 


ing the new method were compared with those obtained using 
the total count method presented in the literature [25]. Ac- 
cording to [23], the second set of counting intervals was se- 
lected for application to the total count method. The first and 
second sets of counting intervals were (1,4), (5,20), (21,40), 
(150,250), (360,560) and (2,5), (6,20), (21,30), (200,300), 
(360,560) min, respectively. The alpha counts from these two 
measurement intervals of the three experiments were input 
into the model of the total count method. Furthermore, a com- 
parison of the results for the activity concentrations of 7”’Rn 
and ??°Rn progeny in air, calculated via the new method and 
total count method, is shown in Fig. 5. 

Fig. 5 illustrates that the activity concentrations of 7”Rn 
and Rn progeny, calculated using the new method, are all 
positive, whereas the activity concentrations of RaB, RaC, 
and ThB, calculated using the total count method with two 
sets of counting intervals, exhibited negative values. The 
occurrence of negative concentration of 7’Rn and 7°Rn 
progeny is not in accordance with reality. The primary rea- 
son may be that the accuracy of calculating comparatively 


lower progeny concentrations in the air using the total count 
method is insufficient, coupled with the statistical fluctua- 
tions in the counts as well as the overlapping effect between 
different spectra during long-term measurements. The total 
count method in [25] utilizes the radioactive decay law of 
alpha particles within five different time intervals (560 min) 
to obtain the concentrations of ?”?Rn and Rn progeny in 
the air. Compared to the total count method, the advantage 
of the new method lies in distinguishing nuclides with over- 
lapping spectra by utilizing the alpha particle counts of non- 
overlapping spectra within consecutive measurement cycles 
(90 min). This enables the distinction of nuclides with dif- 
ferent energy spectra and facilitates nonlinear fitting based on 
the radioactive decay law to determine the concentration of 
each progeny. Moreover, the measurement time is reduced by 
84 

Given the close proximity of the alpha energies emitted by 
RaA and ThC, both around 6 MeV, their counts were not dis- 
tinguished in each experiment, but were counted together. For 
alpha counts with energies of 6 MeV, the quantity of RaA 
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Fig. 5. (Color online) Three experimental results of simultaneous measurements of Rn and “Rn progeny concentrations: (a) first exper- 
iment; (b) second experiment; and (c) third experiment. (The counting interval of the new method was 90 min; the first and second set of 
counting intervals are (1,4), (5,20), (21,40), (150,250), (360,560), and (2,5), (6,20), (21,30), (200,300), (360,560) min) . 


will decay rapidly within a short period because of the much 
shorter half-life of RaA when compared to that of ThC. When 
the total count method with long time intervals is employed, 
the activity concentration of RaA in air may have been under- 
estimated, leading to an overestimation of the activity concen- 
tration of ThC. However, by utilizing the continuous counts of 
the 8.78 MeV ThC’ spectrum within consecutive cycles and 
the determined decay-branching ratio, it is easier to distin- 
guish between the spectral counts of RaA and ThC. Hence, 
the concentration of RaA, calculated using the total count 
method, is lower than that determined by the new method, 
whereas ThC exhibits a higher concentration when compared 
to the calculation by the new method, as shown in Fig. 5. 
Although the fitting quality of some curves in Fig. 4 is not 
optimal, it has minimal effect on the overall accuracy of the 
experiments. Detailed data are listed in Table 2. 


Considering that the relative uncertainty of the detection 
efficiency og /E is + 2.0%, the relative uncertainty of the 
flow rate oQ /Q is +3.0%, and the statistical uncertainty, rel- 
ative uncertainty of radon progeny concentrations provided 
by the new methods are all within +16%, and those of the 
thoron progeny concentrations are within + 13%. These un- 
certainties are significantly lower than those calculated using 
the total count method. Owing to the lower concentrations 
of measured progeny, there are significant statistical fluctua- 


tion errors, but the statistical errors of progeny concentrations 
calculated by the new method are within a reasonable range 
using the 10-min average as the measurement value for one 
cycle. A low measurement uncertainty is more meaningful 
for accurate measurements. 

Given the negative values observed in the concentration of 
?22Rn progeny calculated using the total count method in the 
three experiments, the equilibrium equivalent radon concen- 
tration was not considered. The ?”°Rn decay product concen- 
trations were measured in terms of equivalent concentrations 
(EECra), as provided by the following relations [43]: 


EEC ty = 0.913 (C212p,,) + 0.087 (C214p;) (18) 


where C212p,, C214p; are the activity concentrations in 
Bq/m?. 

Table 3 shows that the Æ EC'rn trends calculated by both 
methods are consistent, with EEC y, in the second exper- 
iment being slightly higher than those in the first and third 
experiments. The lower limit values of the standard errors 
calculated using the total method for the three experiments 
were negative, which is inconsistent with the actual situation. 
However, the values of EEC, calculated using the new 
method are within the range of the standard errors of those 
calculated using the total method, with the relative standard 
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Table 2. Comparison of results from two methods in experiments of simultaneously measuring Rn and **°Rn progeny concentrations (k=1). 


Activity concentration 


Activity concentration 


Activity concentration 


Experiment Progeny (the new method) (the total count method with (the total count method with 
ailai (Bq: m~?) first set of time internals) (Bq - m7?) second set of time internals) (Bq - m~?) 

1 NRaaz 17.234+1.537 0.0090.006 0.002+0.005 
NraB2 3.845+0.364 -37.547424.026 -12.112+33.835 
Nrac2 2.809+0.362 -36.851+24.640 -11.598+35.417 
Nrnp2 2.539+0.199 -1.521+2.850 1.915+2.564 
Nrnc2 1.298+0.124 47.867+27.200 19.166+38.491 

2 Nraa2 23.175+3.495 0.007+0.008 0.008+0.007 
NraB2 7.569+0.491 -22.951+31.027 -47.278+46.204 
Nrac2 6.098+0.384 -23.383431.980 -47.063+48.296 
Nrnp2 4.997+0.386 1.601+3.693 1.298+3.497 
Nrnc2 1.729+0.225 37.305+35.205 63.546+52.534 

3 Nraa2 20.665+3.453 0.010+0.006 0.003+0.005 
NraB2 4.006+0.446 -36.594+25.362 -18.219+36.999 
Nrac2 3.178+0.338 -35.076+25.994 -14.361+38.546 
Nrnp2 3.155+0.161 -1.552+3.007 1.499+2.795 
Nrnce2 1.146+0.088 47.124+28.702 24.997+41.981 


Table 3. Comparison of results of EEC'y, from two methods (k=1). 
EEC rn calculated by 


EEC rp calculated by 


Experiment the total method 
the new method : : 
number (Bq - m~?) with second set of internals 
f (Bq: m~’) 
1 2.431+0.194 3.416+11.615 
2 4.7130.375 6.714+15.851 
3 2.9800. 156 3.5434 12.667 


deviation of EECy, ranging from 5% to 8%. Given the mea- 
surement uncertainties associated with different methods and 
fluctuations in progeny concentrations across various sam- 
pling processes, it can be confirmed that consistency was 
achieved. The experiments were conducted at relatively low 
progeny concentrations of 7”*Rnand ”°Rn. The results ob- 
tained by comparing the two methods demonstrate that the 
proposed method has significant advantages. 


IV. CONCLUSIONS 


The concentrations of ”?Rn and ?°Rn progeny in the en- 
vironment pose a risk of natural radiation exposure to the 
public. A novel mathematical model, utilizing a nonlinear 
fitting method based on the radioactive decay law, is pro- 
posed for the simultaneous measurement of 7?’Rn and 7°Rn 
progeny concentrations in air. The experimental findings in- 
dicated that Rn and ”°Rn progeny concentrations deter- 
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